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ABSTRACT 

Context. The two nuclei of the starburst galaxy Arp220 contain multiple compact radio sources previously identified as radio super- 
novae or supernova remnants. 

Aims. In order to search for an embedded radio AGN, or other possible exotic objects, we have carried out a program of VLBI moni- 
toring at 6 cm over three epochs each separated by four months. 

Methods. Combining the new data with existing data at 6 cm and 18 cm (spanning 4 and 12 years respectively) we are able to charac- 
terise source flux density variability on a range of time-scales. Additionally we analyse the variability of sources in shape and position. 
Results. We detect rapid (< 4 months) variability in three sources (W7, W26 and W29). These sources show possible superluminal 
motion (> Ac) of jet-like features near rapidly varying almost stationary components. These enigmatic sources might be associated 
with an AGN or a highly beamed microquasar (i.e. microblazar). Other hypotheses include that the apparent variability is intrinsic 
and is produced by neutron star powered central components within a supernova remnant, by a sequence of several supernovae within 
super star clusters, or is extrinsic and is produced by Galactic interstellar scintillation of very compact non- varying objects. 
Conclusions. A microquasar/microblazar origin seems to be the best explanation for the nature of the variable sources in Arp220. 
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1. Introduction 

Arp220 is the result of the merger of two galaxies (Norris 1988) 
and has two distinct nuclei separated in projection by ~ 364 pc 
(Scoville et al. 1998). It is the closest ULIRG (Ultra Luminous 
Infra-Red Galaxy) located 77 Mpc away with a far infrared lu- 
minosity Lfir ~ 1.3 x 10 12 L© (Soifer et al. 1987), more typical 
of star-forming galaxies at redshift z = 1. It has a similar star- 
formation density per unit area as redshift z = 6 proto-galaxies 
(Walter et al. 2009). Smith et al. (1998) made the first detection 
at 18 cm of compact radio objects associated with the starburst. 
These objects were detected at shorter wavelengths (13, 6 and 
3.6 cm) by Parra et al. (2007). Both this paper and more recently 
Batejat et al. (2011) concluded that most compact sources are 
either radio supernovae (SNe) or supernova remnants (SNRs). 
However, not all sources were easy to classify. In this paper we 
concentrate on three variable sources detected in the Western 
nucleus of Arp220. In Section 2 we present the observations and 
details of the data processing. In Section 3 we describe our re- 
sults. In Section 4 we discuss several hypotheses for the nature of 
the three sources. Finally, in Section 5 we give our conclusions. 

2. Observations 

We performed deep global VLBI 6 cm observations at three 
epochs (experiments GC031A, GC031B and GC031C) in June 



2008 (epoch 1), October 2008 (epoch 2) and February 2009 
(epoch 3) with the goal of looking for short time- scale variabil- 
ity in sources of all kinds, i.e. due to type Ib/c supernovae or 
AGN candidates. In order to limit bandwidth and time-average 
smearing effects we chose a time resolution of 2 s with frequency 
resolution of 1 MHz. We processed these data in a similar way 
to Batejat et al. (2011). All epochs have been phase-referenced. 
The absolute astrometric accuracy is estimated to be ~ 0.5 mas 
(Batejat et al. 2011). By combining the three epochs (each with 
noise rms ~15 /Jy/beam) we have obtained a very high sensitiv- 
ity image with a noise rms of about 8 yuJy/beam (Figure 1). 

3. Results 

Comparison of our three new 6 cm epochs uncovered three 
rapidly variable sources, W7, W26 and W29 (Batejat et al. 
2010). During 16 years of 18 cm wavelength VLBI monitor- 
ing W7 and W26 have always been present (resp. named SN2 
and SN9 by Rovilos et al. 2005). W29 has been detected dur- 
ing the last 8 years of high sensitivity observations (Lonsdale 
et al. 2006) and the data are consistent with it being present at 
the same flux density before that. All three sources appear to be 
persistent, long lived sources. Those sources, all located in the 
Western nucleus of Arp220, also appear to vary rapidly in posi- 
tion, shape and flux density. Images of these sources at the three 
epochs of project GC031 are shown in Figure 2 together with 
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Fig. 1. Greyscale shows the combined 5 GHz global array GC031A,B,C natural weighted image of the Western nucleus. The noise rms is 
8 /Jy/beam, the image size is 512 x 256 mas. Labeled in white are the variable sources W7, W26 and W29. The green crosses give the 1.3 mm and 
2.6 mm IRAM positions of the hot dust "AGN" feature from Downes & Eckart (2007), the green circle is the position from the 0.8 mm observation 
of Sakamoto et al. (2008). Errors on these positions are estimated to be of order 50 - 100 mas. The black ellipse gives the orientation and size of 
the hot dust feature as fitted by Downes & Eckart (2007) centred at the centroid (red square) of the three millimeter interferometer positions. 



W18, a resolved SNR candidate, as a control source. The light 
curves of those four sources at 18 cm (1994.87 - 2006.43) and 
at 6 cm (2006.02 - 2009.16) are shown in Figure 4.2. Gaussian 
model fitting was applied in the image plane to characterise the 
source structural variations. In some cases somewhat better fits 
could be achieved by fitting for more than two Gaussian com- 
ponents or allowing a component to be significantly extended in 
one direction. In all cases however one or two dominant compo- 
nents were required and their positions and total flux densities 
were robust independent of the starting model. The parameters 
of these components are reported below. 

W7 contains in all epochs a dominant component with al- 
most constant position. In epoch 1 it had a total flux density of 
261 yuJy. In epoch 2 this feature's total flux density decreased to 
174 yuJy while in epoch 3 its total flux density increased again 
to 254 //Jy. In epoch 2 a secondary feature appeared ~1.2 mas 
away to the East (corresponding to - 0.44 pc in Arp220) with 
total flux density 287 //Jy. To obtain the best possible fit to the 
image we require that this secondary component must either be 
resolved or a third component must exist between the main com- 
ponents. Nothing is detected above the noise of the secondary 
component at the same position in epochs 1 or 3. 

W26 contains in epoch 1 a single dominant component with 
total flux density 269 yuJy. In the subsequent two epochs a bright 
component persists at the same position with total flux density 
of respectively 271 fiJy and 442 yuJy. In epoch 2 a possibly re- 
solved secondary feature briefly appears ~1.9 mas to the South- 
West (corresponding to ~ 0.7 pc in Arp220) with total flux den- 
sity 218 yuJy. Nothing is detected above the noise at the same 
position in epochs 1 or 3. 

W29 contains a single component in epoch 1 with total flux 
density 340 //Jy. This component rapidly decreases in flux den- 
sity in subsequent epochs leaving only a weaker feature with to- 
tal flux density 172 //Jy and 140 //Jy respectively, located at ap- 
proximately the same place. In epoch 2 a much brighter feature 
lies ~1.9 mas to the East (corresponding to ~ 0.7 pc in Arp220) 
with nothing visible at this position in the earlier epoch. This 



feature has total flux density 226 /Jy. In epoch 3 this component 
maintains its position and flux density (201 yuJy). 

4. Discussion 

The behaviours of W7 and W26 are similar in the sense that 
they appear to have a persistent main component varying in flux 
density with an additional component appearing only at epoch 2. 
W29 has a single component at epoch 1 which fades away with 
a new component rising and becoming dominant at subsequent 
epochs. If the structural variations are due to components that are 
ejected these results imply apparent motions with velocities > Ac 
for W7 and > 6c for W26 and W29. These velocities are similar 
to the velocity measured in the recently discovered object in M82 
(Muxlow et al. 2010). In Arp220 however the sources seem to 
be persistent while the source in M82 is new. Furthermore, the 
source luminosity of ~ 10 27 erg s" 1 Hz -1 in Arp220 is a factor 
-100 times larger than the luminosity of the source in M82. 

Alternatively, instead of component motion we could be 
observing stationary components that independently rise and 
fade. Flux variations could be intrinsic but also due to extrin- 
sic Galactic scintillation. At our observation frequency, at the 
high Galactic latitude of Arp220 (b = 53°) the critical angu- 
lar scale below which a source shows deep modulation in flux 
density is - 6 yuarcsec and the predicted timescale is of order 
2 hours (Walker 1998). For - 200 /Jy components scintilla- 
tion can plausibly contribute to variability. If the component is 
synchrotron self-absorbed, peaks at 5 GHz and has approximate 
magnetic field and electron energy equipartition, a source diam- 
eter of - 7 yuarcsec is predicted (Chevalier 1998). Even smaller 
source sizes are predicted if the source is relativistically beamed 
(Readhead 1994). 

We note that the hints, found in model-fitting (see Section 3), 
that some sources may require at some epochs either more than 
two compact components along a line or are extended is an argu- 
ment against the observed structural variability being explained 
entirely by scintillation or intrinsic variability of stationary com- 
ponents. We consider however that the quality of our data is 
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not yet good enough to completely rule out these explanations. 
Given this, we consider below physical explanations utilising all 
three variability mechanisms i.e. intrinsic and extrinsic flux vari- 
ability and emission of jet components. 
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Fig. 2. Variability of W7, W26 and W29 over time at 5 GHz. The 
columns show in order the varying sources W7, W26 and W29 and 
finally as a control the non- varying resolved source W18 at epochs 1 to 
3 (from top to bottom). The contours are at 4, 6, 8, ... times the noise 
rms for the variable sources and at 4, 12, 20, ... times the noise rms for 
W18. The noise rms is approximately 15 yt/Jy/beam at all 3 epochs. The 
beam is plotted for each source on the last row. 



4.1. Supernova remnants and supernovae 

Variability within an SNR could be extrinsic and be due to com- 
pact scintillating knots (see introduction to Section 4) in an 
SNR shell. However this would require that a large fraction of 
the flux density of the SNR is contained in very small knots. 
Observations of the M82 SNRs (Fenech et al. 2010) show a 
much smoother flux density distribution over the whole SNR 
shell than required to explain the variability in the sources in 
Arp220. Alternatively variability within SNRs could be intrin- 
sic and originate from plerion central components powered by 
neutron stars (Bietenholz & B artel 2008). 

Another way to explain intrinsic variability of the sources 
in Arp220 would be with the rapid evolution of multiple SNe 
in super star clusters (SSCs). Type Ib/c SNe have radio rise 
times of < 1 month and can reach the luminosity required to 
give 100 //Jy sources at the distance of Arp220 (Chevalier et al. 
2006). Bondi et al. (2012) suggest that source A27 in Arp299 
is such a Type Ib/c SN with rise and decay times < 6 months. 
Assuming super star clusters with masses ~ 10 8 M© equal to the 
most luminous stellar cluster known to date, W3 in NGC 7252 
(Maraston et al. 2004), then Starburst99 (Leitherer & Heckman 
1995) gives, assuming an instantaneous burst, a maximum total 
SN rate of ~ 0.1 yr" 1 , well below the rate of ~ 3 yr" 1 of SN type 
Ib/c needed to explain our observations. More massive SSCs in 
Arp220 are ruled out given that the total dynamic mass of the 
Western nucleus is estimated to be below 10 9 M (Sakamoto 
et al. 1999). 



4.2. Active Galactic Nuclei 

It is possible that at least one of the variable sources is associ- 
ated with a buried radio-loud AGN within Arp220 similar to that 
discovered by Perez-Torres et al. (2010) in Arp299. The AGN in 
Arp299 has a luminosity at 5 GHz of ~ 1.8 x 10 27 erg s" 1 Hz" 1 
comparable to our variable sources. For Arp220 Downes & 
Eckart (2007) list a number of observations from X-ray to radio 
consistent with an AGN in the Western nucleus. Additionally 
they presented dynamical evidence for such a SMBH in the 
Western nucleus of Arp220 by analysing the position-velocity 
diagram from CO(2-l) observations. Downes & Eckart (2007) 
and Sakamoto et al. (2008) have analysed the mm continuum 
emission from hot dust and have argued that it may be too com- 
pact to originate from a starburst and must instead be black hole 
powered. The quoted error bars on the absolute position of this 
dust feature are large; despite this it is interesting that the mean 
of these positions (see Figure 1) lies close to W26 and W29. 
Those two sources are possible AGN candidates whose apparent 
structural variability could be due to the combination of super- 
luminal motion in jet components and scintillation of core com- 
ponents. An obvious problem for the supermassive black-hole 
(SMBH) interpretation is that we have three variable sources. 
Furthermore, if beaming is required to explain the variability an 
even larger parent population is required to explain the presence 
of three detectable SMBH in Arp220. One way to ease this prob- 
lem would be to explain variability within W26 and W29 by the 
scintillation of multiple hotspots in a compact symmetric object 
(CSO) made up of the two sources. The separation between W26 
and W29 is only ~ 7 pc. Most CSOs have sizes around 100 pc 
(Augusto et al. 2006). CSOs of the small size of the system 
W26/W29 would be unusual but not unprecedented. For com- 
pleteness we note that the probability that at least one variable 
source is a background AGN is only ~ 5 x 10" 5 (estimated based 
on the detection of 3 compact radio sources of comparable flux 
densities in the Hubble Deep Field, Garrett et al. 2001). 

4.3. Mlcroblazars 

In this section we consider whether we could explain the vari- 
ability in the three sources assuming that they are beamed mi- 
croquasars (i.e. microblazars). Microquasars are accreting stellar 
mass (-10 M©) black holes. An extragalactic microquasar ori- 
gin has previously been suggested by Muxlow et al. (2010) and 
Joseph et al. (2011) for the superluminal source in M82 while 
Bondi et al. (2012) suggests the same origin for a rapidly vari- 
able source in Arp299. Assuming core-jet sources with T = 5 
observed at an angle to line-of-sight (LOS) of 1 /r, which is the 
typical angle for the brightest selected sources in a randomly ori- 
ented population, the de-beamed intrinsic (rest-frame) radio lu- 
minosity of the Arp220 sources would be ~ 7 x 10 33 erg s" 1 (cal- 
culated by integrating up to the observed frequency and assum- 
ing a flat spectrum). This is still much larger than the mean hard- 
state radio luminosity of typical microquasars (~ 2xl0 30 erg s" 1 , 
Gallo et al. 2012) but these have low accretion rates relative to 
their Eddington accretion rates (< 0.001). If instead we consider 
the extremely luminous Galactic microquasar, GRS 1915+105, 
which accretes at a high fraction of its Eddington rate (Rushton 
et al. 2010) and consider its most luminous radio state, then the 
required rest frame luminosity can be reached. Correcting for its 
observed superluminal motion (Fender et al. 1999) and the re- 
sulting de-boosting of its jet due to it being orientated close to 
the sky plane, a rest frame 6 cm luminosity of ~ 3.4 x 10 33 erg s" 1 
is estimated. Assuming r = 5 and changing its orientation to an 
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Fig. 3. Light curves of W7, W26 and W29 at 18 cm (blue) and 6 cm (red). The SNR candidate W18 is shown as a control source. Error bars are 
plotted at ±l<x. 



angle 1/T to the LOS, we find that GRS 1915+105 would have 
an apparent radio luminosity of about half that of the sources 
in Arp220; for a slightly smaller angle of 0.8/T their luminosi- 
ties are matched. Based on this assumption and under the mi- 
croquasar hypothesis, the double structure observed in the three 
variable Arp220 sources would be explained as a beamed core 
and an approaching beamed component. For a T = 5 source ob- 
served at an angle to LOS 1 /T the apparent superluminal motion 
of the latter is of order 5c, consistent with observations (see in- 
troduction to Section 4). 

A potential problem with the microblazar hypothesis is the 
large apparent size of the double sources in Arp220 (0.44 to 
0.7 pc) and their long fading times (< 4 months) which are 
both much larger than for GRS 1915+105 (0.015 pc and tens 
of days respectively (Fender et al. 1999)). For a GRS 1915+105 
source seen end-on, the discrepancy in scales/timescales is fur- 
ther exacerbated by the effects of geometric projection and rela- 
tivistic time compression. However, not all microquasar jets are 
as compact and rapidly varying as in GRS 1915+105. Hao & 
Zhang (2009) show moving components with r = 3 in the mi- 
croquasar XTE J1550-564 and other sources extending out to 
0.5 pc and persisting for years. If pointed toward us, component 
decay timescales would be relativistically compressed down to 
the order of a few months. Linear scales would still be smaller, 
but of the same order of magnitude, than those observed in the 
Arp220 varying objects. Another possible explanation for du- 
plicating the source sizes observed in Arp220 would be to have 
unresolved beamed GRS 1915+105 like objects within clusters; 
individual sources would vary either intrinsically or because of 
scintillation due to foreground ionised gas in our Galaxy (see 
introduction to Section 4). 

If we assume that only beamed X-ray binaries with jets hav- 
ing an angle to LOS < 1/T can be detected in the radio, then 
W7, W26 and W29 represent only 1% of the parent population 
of unbeamed X-ray binaries and therefore ~ 300 of these are ex- 
pected in Arp220. The total 2-10 keV soft X-ray luminosity in 
Arp220 is L Arp 220 ~ 7 x 10 40 erg s" 1 (Iwasawa et al. 2001, 201 1) 
corresponding to 225xLgrs_1915+ios (Verrecchia et al. 2007) and 
is therefore consistent with this total population size. Finally as- 
suming a star- formation rate of 200 M© yr" 1 over a continuous 
burst lasting > 10 Myr then at near solar metallicity Linden et al. 
(2010) predict 2000 X-ray binaries brighter than GRS 1915+105 
luminosity, more than enough to explain the required size of the 
parent population. 

5. Conclusions 

We have detected three highly variable sources whose struc- 
tures vary on timescales shorter than 4 months with possible 
superluminal motion. It is difficult to conclusively constrain the 
true nature of those three variable sources. However, although 



not all of the properties of the Arp220 variable objects can be 
matched by any single known Galactic microquasar, their lumi- 
nosity and possible superluminal motion are similar to an object 
like GRS 1915+105 seen end-on. New HSA observations look- 
ing for variability on timescales shorter than 30 days are cur- 
rently being reduced; these will hopefully give us more clues as 
to the nature of these objects. 
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